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Abstract. We investigate the role of size quantization in the vacuum force between 
metallic films of nanometric thickness. The force is calculated by the Lifshitz formula 
with the film dielectric tensor derived from the one-electron energies and wavefunctions 
under the assumption of a constant potential inside the film and a uniform distribution 
of the positive ion charge. The results show that quantization effects tend to reduce the 
force with respect to the continuum plasma model. The reduction is more significant 
at low electron densities and for film size of the order of few nanometers and persists 
for separation distances up to 10-^-50 nm. Comparison with previous work indicates 
that the softening of the boundary potential is important in determining the amount 
of the reduction. The calculations are extended to treat Drude intraband absorption. 
It is shown that the inclusion of relaxation time enhances the size quantization effects 
in the force calculations. 



PACS numbers: 
1. Introduction 

Thin metal films of nanometric size can have physical properties different from their 
bulk counterparts due to quantum size effects (QSE) caused by the valence electron 
confinement in the direction normal to the film surfaces [H [2j |3]. Since thin films are 
among the basic components of modern nanodevices an increasing attention has been 
addressed to the quantum vacuum fluctuations forces (van der Waals and Casimir forces) 
between them. Theoretical determinations of these forces have been based mainly on a 
continuum description of the material dielectric properties, neglecting modifications of 
the electronic structures due to the boundaries, i.e. under the assumption that the film 
dielectric properties be the same of the bulk material pfl El El U\ E] • This assumption 
is not valid for small systems with low electron density and/or large surface to volume 
ratio. In such cases a significant size dependence of the basic properties, like the Fermi 
energy and the density of electron states, arises as a consequence of the electron con- 
finement [9]. 
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In a previous paper [TU] we have shown that the inclusion of QSE in the plasma model 
of a free electron metal leads to modifications of the force intensity between nanometric 
size films with respect to the bulk plasma model which range from several to few percent 
depending upon the film electron density and the separation distance. The calculated 
forces show quantum size oscillations and are less intense compared to those determined 
with the bulk plasma model. The calculations were performed using the particle in a 
box model [11] (hereafter indicated as PBM) in which independent electrons are con- 
fined along the direction normal to the film surfaces by hard walls and behave as a two 
dimensional gas parallel to the surfaces. Such model represents a simplified picture of 
the one electron potential in the film. 

In this paper we improve our description of the film dielectric properties along two lines 
of development: first we introduce a finite well model for the one-electron potential 
along the surface normal, second we include intraband absorption by introducing the 
relaxation time in a manner that allows to keep number conservation. As to the first 
issue, we notice that the use of a soft confining potential like a finite well, besides being 
closer to the real shape of the one-electron potential as determined by first principles 
calculations [T2l [TBI [HI EE], allows for a better treatment of QSE for films of small 
thickness, since it does not introduce a priori a distance within which electrons are to 
be confined and leaves the Fermi energy free to oscillate with the film size [9j [16], a 
feature that is not present in the PBM model, where the Fermi energy is kept equal 
to the bulk value. The inclusion of Drude absorption in vacuum force calculations has 
been a subject of considerable debate, mainly in relation with finite temperature cor- 
rections [T71 [HI HH [20l [21j [22]. Very accurate measurements have been reported that 
shows that the simple plasma model of bulk dielectric can give better agreement with 
the experiments than the model with finite relaxation times [23]. Our results refer to 
the T = 0°K case where difficulties do not seem to be present [21]. Our purpose is to 
understand how QSE may affect the force between thin films when intraband absorption 
is included. In this paper we introduce such corrections and we show their importance 
in the calculation of the intensity of the force between metallic films. 



2. Theoretical framework 

We consider two identical metal films of thickness D with plane boundaries separated 
by a distance i (see figure [TJ (a)). D represents the extension, along the z-direction 
normal to the surface, of the positive ion distribution, which is supposed to be uniform 
with the same density of the bulk system. By extending previous results for isotropic 
slabs to the case of films with anisotropic dielectric tensor we can write the expression 
of the force per unit area as [23 [2S1 [23 121] : 

F = ~7T^ kdk du-f(u;) - 7-^ + +^ ~ 7-^ (1) 

2tt 2 J J q LI - Q TM {wjy 1 - Qte{iu) 2 \ 
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Figure 1. (a) Sketch of the two identical interacting films, (b) Infinitely deep quantum 
well with artificial spill-out. (c) Finite quantum well with natural spill-out. 



Q ™ ~ 1 - P 2 TM e-^ M D e Qte ~ i _ p l Ee -2, TE D e ( 2 ) 

_ 7TAf(^) - je xx (u) _ 7tb(w) - 7(w) fQ . 

Ptm — Pte — — [o 

1tm(u) + 7e xx {uj) 7te(w) + 7(0;) 



1te{u) = y A; 2 - ^-e^(^) 7tm(^) = y (77^ ~ ^j) e «H ( 4 ) 




lJ 2 



7 (u;) = WF-- (5) 



c 

Here e xx and e 2 ^ are the diagonal components of the dielectric tensor along the planar 
directions and along the surface normal respectively. We assume e xx = e yy and the off 
diagonal components to be zero. This assumption is consistent with the two dimensional 
gas behavior of the electrons parallel to the surface. The anisotropy of the dielectric 
tensor is a consequence of the finite extension of the film along the z-direction and it is 
the main feature introduced by the size quantization. To calculate the force one needs 
the expression of the dielectric tensor, which has to be derived from the film electronic 
structure. We assume an independent particle model with the one electron potential 
V(z). The electron energies are given by 

£ k„„=^| + ^ (6) 

i.e. they are described by the continuous quantum number ky and by the discrete sub- 
band index n coming from the quantization of the perpendicular wavevector, m being 
the electron mass. The corresponding wavefunctions are given by 

^k ll ,n( r ll^) = -^ eik|rr "^W ( 7 ) 
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here A is the surface area, kn and rn are two dimensional wavevectors parallel to the 
surface, and <p n (z) is supposed to be independently normalized. The functions <p n (z) are 
the solutions of the equation 



I 2m dz 2 



+ V(z)}(j) n (z) = E n <p n (z) (8) 

with the proper boundary conditions. The Fermi energy Ep is obtained through the 
aufbau procedure i.e. by arranging the eigenvalues in ascending numerical order and 
counting until the number of states needed to accommodate all the electrons in the film 
is reached. This procedure leads to a Fermi energy that depends upon the film size 
and is generally different from the bulk value [16]. This can be understood by noting 
that, to ensure charge neutrality, the number of electrons and the number of ions per 
unit area have to be equal. The electron density n(z) can be simply obtained from the 
wavef unctions: 

^) = ^E(^-^)i^)i 2 ( 9 ) 



m=l 



where m is the label of the last occupied state, while the ion density is simply given 
by no = kp B /3ii 2 , where kp B is the bulk Fermi wavevector. By integrating the densities 
along the z axis and imposing that both give the same number of charges per unit area, 
one gets the relation 

^ f2mEF -E n )=n D (10) 



2vr ^ V 



h 2 

m=l 

which, for finite D values, is generally not satisfied if one replaces Ep with its bulk 
counterpart Ep B = h 2 kp B /2m. In the case of the PBM model this equation is not 
satisfied, since one assumes that the Fermi level is the same as in the bulk. To obtain 
charge neutrality one has to impose the additional condition that the electron density be 
confined on a length d larger than D. This artificially introduces the electronic charges 
spill-out but has the consequence that the average electronic density is lower than the 
ion density [3 HUl [2S1 ED]. 

For the purpose of the present study we assume the potential to be a finite well 
V(z) = — Vo inside the film and zero outside. For such finite well model (FWM) the 
energies of the bound states can be written as 

fi 2 h 2 

E v = lL k f + — J&-V (11) 
k lh n 2m 11 2m zn V ' 

where k zn are the quantized transverse wavevectors. They are obtained from the 

equation giving the condition for the existence of bound states in a quantum well [31] : 

W7T 2 i (k zn \ , 
^ = ^-7^ ^ (12) 



'0 



with k = ^2mVo/h. Notice that the first term at the second member is the value of the 
transverse wavevector for an infinite well model of size D and the second term goes to 
zero as Vq goes to infinity. This implies that, for given film size and number of electrons, 
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the Fermi energy referred to the well bottom is higher in the infinite well model. Notice 
that when Vq goes to infinity one does not recover the PBM, since the Fermi energy is 
varied with respect to the bulk value in order to satisfy the charge neutrality condition 
(flQj) . Figured] (b) and (c) illustrate the difference: in the PBM the electronic charge 
density is confined within a distance d = D + A, larger than the size of the ionic charge 
distribution, to allow for the electronic charge spill-out and to ensure global neutrality 
for Ep = Ep B . In the FWM the charge spill-out results naturally from the behaviour 
of the single particle states while the charge neutrality is achieved by varying the Fermi 
energy with respect to the bulk value. Obviously in the limit of infinitely deep well 
(hereafter indicated as IWM) the electronic charge turns out to be entirely localized 
within the length D and the Fermi energy is strongly increased compared to its bulk 
value. 

Once the electron energies and wavefunctions have been obtained, one can calculate the 
dielectric tensor from the expression [TT] : 

up 8ire 2 
u 2 Adm 2 uj^ 



(-') = 1 — - — — E E /(£k||,n)(£k||,n - £k'n') X 



k||,nk'| ,n' 
\{lpk [{ ,n\p a \lpic' r n'}\ 2 
X (£ k|| ,n-£ ki| ,n0 2 -^ 2 (13) 

here a = x, y, z labels the cartesian component of the tensor, p a indicates the component 
of the electron linear momentum, up = Qpn(D) /uq is the plasma frequency of the 
quantized electron gas (f2p = ^ Aire 2 No /m is the free electron plasma frequency and 
n(D) is the average electron density of the film) and /(-Ek^n) is the occupation factor 
of the (k||,n) state. In the IWM and the FWM u P = fi P , while in PBM n(D) is 
smaller than n , since the electronic charge is distributed over a larger distance than 
the ionic charge (see Figfj] (b)). The off-diagonal component are equal to zero. This 
expression differs from the plasma model dielectric function adopted in previous studies 
in that: (i) it has a tensor character with e xx = e yy ^ e 2Z , (ii) through the double 
sum in the second member it accounts for transitions between lateral sub-bands, whose 
probability amplitude is expressed by the momentum matrix element between the one 
electron wavefunctions (J7j). It can be easily shown that these transitions do not affect 
the lateral components of the dielectric tensor, which are given by the simple expression 
of the plasma dielectric function 

U) 2 

e^H = e yy (oj) = 1 - -| (14) 
because the momentum matrix element for x and y component vanishes. 



3. Results for finite well potential 

In this section we present the results of calculations for finite well potentials. We take 
Al, Ag and Cs corresponding to a radius per electron in Bohr units r s /a equal to 
2.07, 3.02 and 5.62 respectively, in order to illustrate QSE at different densities and 




Figure 2. Fermi energy normalized to its bulk value for Al, Ag and Cs using the 
IWM (dashed line) and the FWM (continuous line). The bulk Fermi energies and the 
work functions have been taken from ref. |32j . 



potential depths. The value of the potential depth is obtained by summing the metal 
work function W with the calculated Fermi energy. Figure [2] reports the calculated 
Fermi energies as a function of the product between the Fermi wavevector and the film 
thickness. This allows to better point out the oscillations and the cusps arising from the 
crossing of the Fermi energy by new subbands upon varying the film size. The figure 
shows the V value appropriate to the bulk and to the large D limit. We give in the 
same figure the results obtained by assuming an infinitely deep potential (IWM). The 
comparison allows to illustrate the effects of the potential softening. In agreement with 
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previously published results [TS1 ESI E3] we find that: 

• E F is systematically larger than Ep B and goes to the bulk value as D goes to 
infinity. The difference is more pronounced and the bulk limit is achieved at larger 
size in the low density systems, as it is clearly shown by the comparison between 
Cs and Al curves; 

• As expected from the discussion of the previous section, the softening of the 
potential leads to less pronounced deviations from the bulk values. Because of 
the stronger electron confinement, the IWM has a larger Fermi energy that the 
FWM; 

• The cusps correspond to integer values of half of the Fermi wavelength in the IWM 
case. This feature is only approximately verified for the FWM. 

As a consequence of size quantization, the ^-component of the dielectric tensor is 
expected to go to a finite value e 2Z (0) as the frequency goes to zero. This value increases 
proportionally to D 2 in the large size limit [TTj . In Figure [3] we plot the quantity 
e zz (0)/D 2 for the three cases under study. Again the results show the cusps due to 
the filling of new subbands as D increases [10]. The asymptotic limit is obtained for 
k F D /it of the order of 5 -h 6 in the three cases. Significant differences appear in the large 
D behaviour when the FWM results are compared with those from the IWM model: 
e zz (0) is larger for finite wells of small size, while it is smaller at high D values. The 
convergence to the asymptotic limit is considerably slower for the infinite well, specially 
in the low density metals. This behaviour reflects the differences in the distribution 
of the eigenvalues of equation ([8]). For the infinite well there are infinite bound states 
whose energy scales like n 2 , see equation (fT2|) . and the separation between two successive 
levels increases linearly with n. Such behaviour is not present in the FWM, for which 
equation ([8]) has a finite number of eigenvalues corresponding to bound states and a 
continuum spectrum at positive energies. The PBM curve takes values closer to the 
FWM than to the IWM. This is primarily a consequence of the fact that PBM allows 
for electron charge spill-out, while in IWM the electron distribution is confined within 
D i.e. it has the same size of the positive charge. 

In the following we show the changes in the force caused by the size quantization with 
respect to the results obtained by using the isotropic continuum plasma model, where 
the ^-component of the film dielectric tensor is equal to the planar components i.e. 
is given by equation (fl4l) . We use the symbol Fq to indicate the force per unit area 
calculated for the quantized film, while Fp is the force per unit area calculated in the 
isotropic plasma model. To better illustrate the results, in figures H] and [5] we plot the 
quantity 

F P - F Q 

5 P = - (15) 

fp 

as a function of the separation distance I for films of 1 and 5 nm thickness respectively. 
In each figure we display the results for the three cases under study and we compare 
the finite well with the IWM at the same density. This allows us to point out the 
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Figure 3. Static value of the zz component of the dielectric tensor for Al, Ag and 
Cs using the IWM (dashed lines), the PBM (dotted lines) and the FWM (continuous 
lines). 



modifications caused by the potential softening. We also show the curves appropriate 
to the PBM. In agreement with our previous findings we observe that 

• QSE tend to reduce the intensity of the force; 

• the reduction is more significant at low density (Cs) than at high density (Al); 

• it may be considerably higher than 10% for 1 nm thickness and reduce to few per 
cent at 5 nm; 

• it can be appreciable over a distance interval up to 10 -j- 50 nm. 
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Figure 4. Relative percentual difference for the force between two identical films of 
thickness D = 1 nm as a function of the films separation £, for Al, Ag and Cs. Dashed 
lines have been obtained using the PBM, dotted lines have been obtained using the 
IWM and the continuous lines represents the finite well results. 

The most important conclusion that can be drawn from the figures is that the potential 
shape is important and can lead to a substantial modifications of the quantum size effects 
both at small and at large distances. The models which confine the electronic charge 
tend to overestimate the force reduction induced by size quantization. The curves for 
the PBM show large force reduction (greater Sp values) over a wide interval of distances. 
On passing to the IWM case one notice that the removal of the constriction that the 
Fermi energy be equal to the bulk value, still keeping an infinitely deep potential, leads 
to smaller Sp values and to a more rapid decay of the curves at large distances. Reducing 
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Figure 5. Relative percentual difference for the force between two identical films of 
thickness D = 5 nm as a function of the films separation £, for Al, Ag and Cs. Dashed 
lines have been obtained using the PBM, dotted lines have been obtained using the 
IWM and the continuous lines represents the finite well results. 



the well depth to finite values has a similar effect: it causes a general decrease of Sp 
and a narrowing of the distance interval over which QSE are appreciable. This also 
implies that any increase of the confining potential depth at fixed ion density leads to 
higher 5p values and to more significant QSE. The large values taken by 5p in the PBM 
case do not arise from the charge confinement only, since, as pointed out before, the 
constraint on the charge distribution is weaker that in the IWM. To a large extent they 
are a consequence of the plasma frequency normalization caused by the decrease in the 
average electron charge density that it is necessary in order to achieve global charge 
neutrality [XQl [29J EQ]- In the isotropic plasma model one takes up = Qp. In the PBM 
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this value is obtained only at large film thickness. Neglecting this normalization i.e. 
taking the free electron plasma frequency in the parallel components of the dielectric 
tensor (but not in e zz ) would lead to Sp values closer to the well potential models. 

4. Intraband absorption effects 

To introduce intraband absorption we have to modify the dielectric tensor in a way that 
allows to include relaxation time effects in the parallel components and to recover the 
Drude behaviour in the large D limit. This cannot be done by simply introducing an 
imaginary part of the frequency u, since this violates the continuity equation locally 
[341 135] . The appropriate recipe is to replace into equation (1T31) uj 2 with uj{uj + z27t/t), 
where r is the relaxation time. For the parallel components this leads to the Drude 
dielectric function 

= 1 + , (16) 

io(u + 27) 

where 7 = 2-k/t The results of the calculation of the force per unit area with the 
modified dielectric tensor are displayed in Figures [6] and [7] for 1 nm and 5 nm films. We 
report the quantity 

Fd{i) 

where Fj) is the force calculated using the bulk Drude model with a given relaxation 
time and Fqd in the one obtained by the calculation with the same relaxation time 
and with size quantization included. For r — > 00 we recover the plasma model so that 
8p = <5d(0). The calculations have been performed assuming a finite potential well and 
for two values of the relaxation frequency 7. We have taken values that approximately 
correspond to those reported for the metals under consideration [32]. The figures show 
a comparison with the curves obtained by the plasma model with QSE. It is clear from 
these results that the main effect of the inclusion of intraband absorption is to increase 
5 i.e. to increase the difference with respect to the calculations with the bulk dielectric 
function. The smaller is the relaxation time the larger is the reduction of the force. The 
effect is qualitatively the same in the three metals under study, but it depends upon the 
well depth and the film size. The influence of the shape of the potential is illustrated 
in FigJH] showing the curves of Ag films at a given relaxation frequency for the various 
models. Again it should be noticed that the FWM gives the lower 5d values. The PBM 
results show large 5d values over a very wide interval of distances. To a large extent this 
behaviour has to be imputed to the renormalization of the plasma frequency. Figure [9] 
shows typical curves of 5n as a function of the film thickness for different values of the 
relaxation frequencies at a given separation distance of 5 nm. As expected 5d decreases 
with D, but the slope at large thicknesses (D of the order of 10 4- 50 nm) depends 
significantly upon the relaxation frequency. 




Figure 6. Relative percentual difference for the force between two identical films 
of thickness D = 1 inn as a function of the films separation £, for Al, Ag and Cs. 
Continuous lines represent the results for 7 = 0. Dotted lines have been obtained 
using 7 = 5 x 10 13 rad/s for Ag and Cs and 7 = 10 14 rad/s for Al. Dashed lines have 
been obtained using 7 = 10 14 rad/s for Ag and Cs and 7 = 10 15 rad/s for Al. 



5. Conclusions and lines of development 

We have presented a rather complete set of theoretical results illustrating the possible 
role of size quantization effects in the electromagnetic vacuum force between very thin 
films and showing how the determination of these effects depends upon the description 
of the film electronic structure and upon the inclusion of intraband effects. However 
we want to point out that our analysis is still a mere indication of the corrections 
to the simple picture that assumes the same dielectric function for films and bulk 
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Figure 7. Relative percentual difference for the force between two identical films 
of thickness D = 5 nm as a function of the films separation £, for Al, Ag and Cs. 
Continuous lines represent the results for 7 = 0. Dotted lines have been obtained 
using 7 = 5 x 10 13 rad/s for Ag and Cs and 7 = 10 14 rad/s for Al. Dashed lines have 
been obtained using 7 = 10 14 rad/s for Ag and Cs and 7 = 10 15 rad/s for Al. 



solids. Although the basic features of size quantization (confinement of the electronic 
charge, anisotropy of the dielectric tensor, presence of inter-subbands transitions in the 
dielectric function) are already present in the models we have studied, there is room for 
substantial improvements before an accurate comparison with experimental data, like 
those obtained by Lisanti et al. for Pd films [36], can be done. A more detailed treatment 
should include (i) band structure effects, (ii) non-locality of the dielectric response and 
(iii) non-local treatment of the reflectivity. A self consistent first principles calculation 
of the inverse dielectric matrix for a slab of the appropriate size and symmetry, from 
which a macroscopic dielectric function can be derived with band structure and non- 
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Figure 8. Relative percentual difference for the force between two identical Ag films of 
thickness D = 5 nm as a function of the films separation £, 7 = 10 14 rad/s. Continuous 
lines represent the FWM, dotted line the IWM and dashed line the PBM. 
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Figure 9. Relative percentual difference for the force between two identical Ag films 
separated by a distance t — 5 nm as a function of the films thickness D. Continuous 
line represents the result for 7 = whereas the dashed line has been obtained with 
7 = 10 14 rad/s. 

local effects included, could provide the appropriate treatment of the first two issues 
[37J [38]. Corrections to the Fresnel optics, along the lines indicated by several authors 
(39J HOI HI] , can lead to important modifications of the reflectivity even in the case of a 
free electron gas film. We are currently investigating these matters and the results will 
be presented elsewhere. Still when comparing theory with experiments for thin films one 
should consider the fact that measured relaxation frequency turn out to depend upon 
the film size and morphology [121 S3J, HH HS1 H6] . This is largely due to the so called 
classical size effects arising from the scattering of the electrons at the film boundaries. 
In view of the sensitivity of QSE to the value of the relaxation time, comparison with 
experimental data may be possible only if a realistic estimate of the modifications in 
the relaxation times caused by surface scattering is available. 
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